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We have studied damping in poly crystalline Al nanomechanical resonators by measuring the 
temperature dependence of their resonance frequency and quality factor over a temperature range 
of 0.1 - 4 K. Two regimes are clearly distinguished with a crossover temperature of 1 K. Below 
1 K we observe a logarithmic temperature dependence of the frequency and linear dependence of 
damping that cannot be explained by the existing standard models. We attribute these phenomena 
to the effect of the two-level systems characterized by the unexpectedly long (at least two orders 
of magnitude longer) relaxation times and discuss possible microscopic models for such systems. 
We conclude that the dynamics of the two-level systems is dominated by their interaction with 
one-dimensional phonon modes of the resonators. 

PACS numbers: 85.85.+j, 62.25.-g, 81.05.Bx, 62.40.-H 
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Nanoelectromechanical systems have recently attained 
a lot of interest due to a variety of promising applica- 
tions such as ultra-sensitive mass measurement |l| and 
single-spin detection Q as well as due to their suitabil- 
ity for doing quantum- limited measurements ['ij . For 
this, nanomechanical resonators with high resonance fre- 
quencies in combination with high quality factors are re- 
quired. So far, nanomechanical resonators mostly made 
of single crystal semiconductor and dielectric materials 
have been studied, and frequencies above 1 GHz with 
quality factors of about 500 have been achieved @ . Low- 
temperature studies on such resonators reveal glass-like 
behavior 7. 8] and suggest that this may be due to the 
interaction of the flexural beam modes with two-level- 
systems (TLS). Metallic beams can be fabricated using 
technologies known for metallic nanoelectronic devices, 
which therefore allows easy integration of mechanical de- 
grees of freedom into such devices. Although fabrication 
methods for metallic beams of nanoscale size have been 
developed recently 0, [l(| , very little is known about loss 
mechanisms in such beams. In this respect, understand- 
ing damping in nanoscale mechanical resonators is of pri- 
mary importance. In this work we demonstrate that the 
quality factor of Al doubly clamped beams at mK tem- 
peratures is affected by two-level systems. The observed 
temperature dependences of damping and resonance fre- 
quency are characteristic for amorphous materials with 
TLS relaxation due to one dimensional phonons. Such 
a dependence cannot be interpreted within the standard 
TLS model ll|, |12[ , and we propose an alternative ex- 
planation. The apparent similarity between amorphous 
insulators and polycrystalline metals and their dramatic 
difference from amorphous metals implies a different na- 
ture of TLS in metallic beams. 

Doubly clamped beams are fabricated on an oxidized 
silicon substrate using a trilayer resist structure with an 



extra sacrificial calixarene layer (lol . 13 1 . The beam mask 
in Ge layer is defined by electron-beam lithography and 
reactive ion etching. After metallization and lift-off pro- 
cess, the beams are suspended by removing the underly- 
ing calixarene layer in an oxygen plasma. The Al poly- 
crystalline beams, with length I, width w and thickness t 
(see left inset of Fig. [T]), are connected to the central line 
of prefabricated gold-patterned on-chip coplanar waveg- 
uides which are ribbon-bonded to coaxial lines. The sam- 
ples are mounted in vacuum space inside the bore of a su- 
perconducting solenoid providing a transverse magnetic 
field of up to 5 T. All measurements were done in a dilu- 
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FIG. 1: (color online) Amplitude of the induced electromotive 
force of the 5 fim beam at B = 0.5 T for different temperatures 
(symbols). The solid lines are Lorentzian fits. Right inset 
shows damping 1/Q of the same beam as a function of the 
magnetic field (black squares). The independence expected 
from magnetomotive damping is shown as a fit (red line) to 
the experimental data. Left inset presents a layout of the 
beams studied. 
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tion refrigerator with a base temperature down to about 
50 mK. 

We characterize our beams usin g a conventional mag- 
netomotive measurement scheme [14j. The RF signal 
from the output of the network analyzer is fed into the 
coaxial line at the top of the cryostat and delivered to 
the chip through a 20 dB attenuator at the 4 K stage. 
AC current flowing through the beam in a perpendicular 
external magnetic field B actuates the device due to the 
Lorentz force. On resonance the beam dissipates energy 
producing a dip in the transmitted signal which is deliv- 
ered via a second coaxial line to the room-temperature 
preamplifier and then to the input of the network ana- 
lyzer. The power applied to the beams is low enough to 
keep them in the linear regime, as seen from the reso- 
nances depicted in Fig. [T]for the 5-/im long beam. Elec- 
tromotive force induced on the beam is estimated from 
the measured transmission as Vemf = V PZ, where P is 
the power difference in the transmitted signal on and off 
resonance and Z = 50 f2 is the impedance of the high- 
frequency line. The resonance frequency /o and damping 
1/Q are extracted from these traces by fitting them with 
aLorentzian. Damping 1/Q — A///o is defined as a full 
width at half maximum divided by the center frequency 
of the Lorentzian fit to the transmitted power. Damp- 
ing dependence on the applied magnetic field is shown in 
the right inset of Fig. [TJ Magnetomotive damping is pro- 
portional to the motional impedance that scales as B 2 , 
therefore it dominates in high magnetic fields and pro- 
duces an overall parabolic dependence. This is confirmed 
by the parabolic fit to the experimental data between 1 T 
and 5 T shown by the red line and extrapolated to lower 
fields. Below about 1 T, damping saturates and deviates 
from the parabolic dependence, which may be attributed 
to the contribution from the field dependent losses in the 
external circuit. 

Figure [2] shows the relative change of the resonance 
frequency <5///g nax , for four beams of different dimen- 
sions, as a function of temperature on a logarithmic scale. 
For each trace, the frequency change is normalized to its 
maximum resonance frequency /™ ax . All the beams ex- 
hibit qualitatively the same overall temperature depen- 
dence: the resonance frequency reaches a maximum at 
about 1.5 K and decreases at higher and lower tempera- 
tures. Moreover, below this temperature, the frequency 
decreases logarithmically down to the lowest temperature 
that could be reached in our measurements, for all four 
beams, although their slopes are different. It should be 
noted here that the beam resonance frequencies change 
with temperature due to tension created due to the dif- 
ference in coefficients of thermal expansion of Al and Si, 
however, this dependence is monotonous in the tempera- 
ture range from 300 K down to about 20 K. Below 20 K 
the thermal expansion coefficients saturate and hence no 
frequency change is expected. Therefore, the character- 
istic frequency dependence at low temperature must be 
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FIG. 2: (color online) Relative shift of the resonance fre- 
quency as a function of the temperature for different res- 
onators (symbols). For each device the frequency shift is 
normalized to the maximum resonance frequency /™". The 
dotted lines are logarithmic fits to the measured data used to 
determine the values for C in Tabled] 



attributed to a different mechanism as discussed below. 

Figure [3] shows damping as a function of temperature 
for the same range as in Fig. [3J but with a linear tem- 
perature scale. Again, all four beams show qualitatively 
the same temperature dependence. Damping increases 
linearly with temperature up to about 1 K. Above this 
temperature it continues to increase but with a signifi- 
cantly smaller slope. 

Now we turn to the analysis of possible loss mecha- 
nisms in Al beams. Damping 1/Q can be divided into 
internal and external as 1/Q = 1/Qmt + 1/Qoxt- 
For the 5-/zm long beam, the observed parabolic depen- 
dence of damping on magnetic field (see inset of Fig. [1]) 
is related to external losses 1/Qcxt in the measurement 
circuit, called magnetomotive damping jljj . Since we 
are interested in the intrinsic or material-dependent loss 
mechanisms 1/Qint, we subtract the effect of magneto- 
motive damping from the measured values of 1/Q for all 
beams. In addition, since magnetomotive damping scales 
as oc all measurements of the temperature dependence 
for the 5-/zm long beam have been performed at 0.5 T. 

Several intrinsic mechanisms discussed in the litera- 
ture can contribute to the dissipation in nanomechanical 
resonators at low temperatures. Thermoelastic damp- 
ing has been shown to be negligible at high frequencies 
and low temperatures and therefore does not play any 
role in the measurements reported here [ToT ]. In metals, 
electron-phonon scattering may also contribute to damp- 
ing [17[ , but its contribution in polycrystalline metals can 
usually be neglected due to the short mean free path of 
electrons [3]. 

Clamping loss is another contribution to damping of 
the beams. It depends on the beam aspect ratio as (w/l) 3 
for in-plane displacement [l9| . therefore this loss mecha- 
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FIG. 3: (color online) Damping \/Q ln ± as a function of tem- 
perature for four resonators of different lengths (symbols). 
Damping increases linearly up to 0.7 — 1.5 K followed by a 
weaker temperature dependence above this temperature. The 
dotted lines are linear fits to the measured data used to de- 
termine the values for vpU ' . 



TABLE I: Physical parameters of Al resonators with different 
dimensions. The beam thickness is 0.2 /im except for the 
narrow beams (w = 0.06 /im) whose thickness is 0.1 /jm. The 
resonance frequency fo and the quality factor Q are extracted 
at T = 100 mK. 



1 X 


w [/im 2 ] 


5x0.1 


3x0.1 


2x0.1 


1x0.06 


0.6x0.06 


fo 


[MHz] 


40.5 


110 


187 


360 


770 




Q 


120000 


41000 


16000 


30000 


2400 




c 


3xl0 -5 


4xl0~ 5 


8xl0~ 5 


17xl0~ 5 


X 




v F U 


0.05 


0.09 


0.08 


0.09 


X 



nism becomes significant for the beams with a high res- 
onance frequency, i.e., small aspect ratio, as can be seen 
in Table HI However, this damping mechanism is not ex- 
pected to show any temperature dependence and there- 
fore contributes to the data shown in Fig. [3] only as a 
constant background without affecting the slope. Ide- 
ally, the temperature dependence of 1/Qint should be 
shifted up for shorter beams provided only their length 
is changed. Such a tendency is observed in Fig. [3] for 
the three beams with the same width and thickness (see 
data for beams with f = 40.5, 110 and 187 MHz). Ta- 
ble I shows that with the fabrication method utilized, a 
resonance frequency as high as 770 MHz together with 
a quality factor of 2400 have already been achieved. By 
optimizing the beam dimensions, it should also be possi- 
ble to push the resonance frequency up to above 1 GHz 
while preserving a reasonably high quality factor. 

Qualitatively, the temperature dependence of the 
resonance frequency and damping in our experiments 
looks similar to those observed in amorphous insulators 
where they are attributed to two- level systems ll|. In 
these materials the temperature dependence of both the 



sound velocity and damping display two regimes with 
a crossover temperature around T* ~ IK. The low- 
temperature regime is characterized by a logarithmic 
increase of the sound velocity with temperature while 
the quality factor changes as a power law with temper- 
ature. The crossover between the two regimes occurs 
when the TLS relaxation rate T becomes approximately 
equal to the measurement frequency uj — 2wf, which 
is close to the beam resonance frequency in our experi- 
ments. One expects that in amorphous metals the inter- 
action of the TLS with electrons is much stronger leading 
to a much larger value of T, so that the high-frequency 
(low-temperature) regime is very difficult to observe, in 



agreement with the existing data [20|, |2l| . The situation 
in polycrystalline metals, such as the samples studied in 
this work, is less clear. In particular, the naive expec- 
tation that the defects in these materials should behave 
similarly to the ones in amorphous metals was not con- 
firmed experimentally. Instead, surprisinglvjthey behave 
more like defects in amorphous insulators [22]. Our data 
also confirm this fact and imply that the interaction of 
the defects with conduction electrons is very small, in 
agreement with Ref. [23j. 

For all mechanisms one expects that the sound velocity 
changes logarithmically in the low temperature regime: 



6 -f = cm(L 

fo \T 



(1) 



where C is a dimensionless parameter that characterizes 
the TLS interaction with sound waves and with each 
other: C — vt'J 2 '/E, where vt is the TLS density of 
states, 7 is the interaction constant and E is the Young's 
modulus. Remarkably, the value of C is known to be al- 
most universal for all amorphous materials, C ~ 10~ 3 - 
10 -4 (24|. Our measurements on polycrystalline Al sam- 
ples give values of C » 0.4-1.7xl0~ 4 , consistent with 
this phenomenology. 

In contrast, temperature dependence of damping is 
very sensitive to the physics of TLS because it directly 
probes the TLS relaxation rate: 
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The observed linear temperature dependence of damp- 
ing below 1 K implies that all data below this tempera- 
ture correspond to the high-frequency regime u) > T(T) 
and that T oc T. For this to be true, even the lowest 
resonance frequency fo ~ 40 MHz must be larger than 
the relaxation rate at T ~ 1 K. Note that the tempera- 
ture dependent part of damping l/Qi n t(T)~ 1/Qint (0) ^ 
0.5 x 10~ 4 is consistent with the regime lu > T(T). 

The conclusion oj > T(T) is difficult to reconcile with 
the electron mechanism of TLS relaxation because for 
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this mechanism HT = 2ir(vpU) 2 k,BT ', where vp is the 
electron density of states and U is their interaction with 
TLS. Assuming that this condition holds for f = 40 MHz 
and T = 1 K one would conclude that vpU < 0.01, which 
is much smaller than one expects and observes for a con- 
ventional TLS in a metal, v F U = 0.1-1 [Til]. The dif- 
ferent temperature dependence of damping studied here 
from that reported for bulk polycrystalline Al samples 
also points towards the phonon dominated relaxation. 

This conclusion, however, is in a perfect agreement 
with the phonon mechanism of TLS relaxation when one 
takes into account the fact that for temperatures T < 1 K 
the phonon wavelength in Al is A > 0.25 //m, so that 
at these temperatures the studied beams are essentially 
one dimensional structures. The linear phonon spectrum 
;25| implies a constant density of states in a ID system, 
and therefore the TLS relaxation rate is expected (cf. 



Ill IJJ, l26j) to be hT ~ (a 2 /wt)k B T, where a is the 
lattice constant. The ratio (a 2 /wt) ~ 10 -5 naturally 
leads to the right order of magnitude of the relaxation 
rate and its linear temperature dependence. Unlike the 
phonon properties, however, the reduced dimensionality 
of our resonators does not affect their electron properties 
or interaction with TLS. 

The apparent absence of the TLS-electron interaction 
in polycrystalline materials points to a different origin of 
the TLS in such materials. In amorphous insulators and 
metals, the TLS are likely to be single atoms that tunnel 
atomic distances between two positions. In cry stals, a 
more likely origin are kinks on dislocations [28(. These 
kinks are very smooth objects in soft metals (such as Au 
and Al) due to a small value of the Peierls barrier and 
thus may interact very weakly with the electrons [29| . 

The conclusions reached above are based on the lin- 
ear temperature dependence of 1/ Q- m t(T) and constant C 
expected for conventional TLS. These assumptions must 
be reexamined for TLS originating from smooth kinks on 
dislocations. For instance, interaction between the kinks 
in strained samples may not be so small as between TLS 
because the external strain creates kinks until the in- 
teraction becomes sufficient to balance the strain. Such 
interaction may suppress the density of states at low ener- 
gies similar to spin glass physics 30]. However, a typical 
model predicting suppression of the density of states and 
a strong linear T-dependence damping would also give a 
sound velocity change different from Eq. (p}. 

Observation of a single TLS and study of its dynamics 
in polycrystalline materials would be the most direct way 
to identify mechanisms discussed above. Alternatively, 
one can probe the interaction with electrons by studying 
the effect of superconductivity on the TLS relaxation rate 
and l/Q int (T). 

In conclusion we have reported measurements of the 
temperature dependence of damping and the resonance 
frequency of the fundamental mode of doubly clamped 
metallic nanomechanical resonators. Our data indicate 



that these are dominated by unconventional TLS with a 
long relaxation time which can be associated with dislo- 
cation kinks. 
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